Evidence for Intra-Unit-Cell magnetic order in Bi2Sr2CaCu208+5 
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Polarized elastic neutron scattering measurements have been performed in the bilayer copper oxide 
system Bi2Sr2CaCu2 08+5, providing evidence for an intra unit cell (lUC) magnetic order inside the 
pseudo-gap state. That shows time reversal symmetry breaking in that state as already reported 
in Bi2Sr2CaCu2 0g+i through dichroism in circularly polarized photoemission experiments. The 
magnetic order displays the same characteristic features as the one previously reported for monolayer 
HgBa2Cu04+4 and bilayer YBa2Cu306+x, demonstrating that this genuine phase is ubiquitous of 
the pseudo-gap of high temperature copper oxide materials. 

PACS numbers: PACS numbers: 74.25.Ha 74.72.Bk, 25.40.Fq 



The existence of the mysterious pseudo-gap (PG) state 
in the phase diagram of copper oxide superconductor and 
its interplay with unconventional d-wave superconduc- 
tivity has been a long standing issue for more than a 
decade. There is now a growing number of experimen- 
tal indications that the pseudo-gap phase actually corre- 
sponds to a symmetry breakingstate [Il-Q . In his theory 
for cuprates, C. M. Varma HjIt!] proposes that PG is a 
new state of matter associated with the spontaneous ap- 
pearance of circulating current (CC) loops within Cu02 
unit cell. This Intra-Unit-Cell (lUC) order breaks time 
reversal symmetry, but preserves lattice translation in- 
variance. While from a theoretical point of view the ex- 
istence of a CC-loop order and the ability of such a q=0 
instability to produce a gap in the charge excitation spec- 
trum are still highly controversial |8l4lOj. several exper- 
imental observations provide strong encouragement for 
models based on CC-loop order in copper oxide materi- 
als. Polarized elastic neutron scatte ring studies carried 
in bilayer YBa2Cu306-i-x (Y123) Ill4l4l| and monolayer 
HgBa2Cu04+5 (Hgl201) [ij, [lef have reported experi- 
mental evidence of a long range 3D magnetic order hid- 
den in the PG state. In La2-xSrxCu04 (Lal24), a sim- 
ilar magnetic order has also been observed but in 
this system it remains 2D and short ranged. This novel 
magnetic state preserves the lattice translation invariant, 
but, at variance with ferromagnets, does not give rise to 
a uniform magnetization [2]. These observations imply 
the existence of an lUC (antiferro-)magnetic order, whose 
symmetry is consistent with the so-called CC-6'ii phase 
proposed by C. M. Varma [1, |3|- Within that model, 
neutron diffraction measures the distribution of static 
magnetic fields generated by the CC-loops. The possi- 
ble detection of these magnetic fields by local probes is 
still debated [l^, Esl , but should help to get a deeper un- 
derstanding of the instrinsic nature of the lUC magnetic 
order. In Y123, the observation of anomalies in the sec- 
ond derivative of the magnetization [20i| when the lUC 
magnetic order settles in and the breaking of time re- 



versal symmetry below the ordering temperature in Kerr 
effect measurements ^ can also be understood within 
the framework of the CC-loop model [21,, 2^, putting 
additional symmetry constraints 23j . In addition, obser- 
vation of CC-loops in Cu02 unit cell has been recently 
reported in non superconducting CuO materials using 
resonant X-ray diffraction 24| . 

Here, we report a polarized neutron diffraction study 
of the bilayer cuprate, Bi2Sr2CaCu08+5 (Bi2212). Up 
to now, the most accurate information concerning the 
electronic properties in the PG state have been provided 
through tunneling spectroscopy [2^ and angle resolved 
photo-emission (ARPES) [l^ carried out mostly in this 
system. Our study shows that the long range lUC mag- 
netic order found in Hgl201 and Y123 is also present 
in Bi2212. This order develops at a temperature T^ag 
close to the PG temperature T* reported by various tech- 
niques. The observation of an lUC magnetic order in 
Bi2212 confirms that time reversal symmetry is broken 
in the PG state, as first suggested by the circular dichro- 
ism in ARPES [H, 0, E^l- Likewise, our polarized neutron 
diffraction study in Bi2212 and the analysis of the STM 
images Q in the same system suggest that lUC order is 
likely involved in the PG physics. 

Measurements are performed on underdoped (UD) 
Bi2Sr2CaCu208+<5 (Bi2212) single crystals. The synthe- 
sis is carried out using the travelling solvent floating zone 
technique (TSFZ) in air [2^. Three large crystals have 
been extracted from the as-grown rod. The composition 
homogeneity and the bulk crystal quality is provided by 
EDX and neutron diffraction studies. The as-grown sin- 
gles crystals are weakly overdoped (OD) with a nominal 
superconducting critical temperature Tc of 87 K. The un- 
derdoping of the samples is then achieved using a post- 
annealing treatment of 300 hours under reduced oxygen 
atmosphere P(O2)=0.05 atm at 450°C, yielding an aver- 
age onset Tc of 85 K [28|. 

For the neutron diffraction measurements, the co- 
aligned single crystals are attached on the cold head of a 
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FIG. 1: Temperature dependence of the intrinsic Bragg scat- 
tering at Q=(1,0,1) in the SF channel (full squares) and of 
polarization leakage (opened circles) corresponding to a T- 
independent bare flipping ratio FR°= 52.6. 



4K-closed cycle refrigerator and aligned in the [100]/[001] 
scattering plane, so that tranferred wave vectors Q of the 
form (H,0,L) are accessible. Q is given in reduced lattice 
units (^7 177^), using tetragonal notations a ~ b = 
3.82 A and c=30.87 A. Polarized neutron diffraction 
measurements are performed on triple-axis spectrometer 
4F1 at reactor Orphee in Saclay (France). The polarized 
neutron scattering set-up is similar to the one used in 
previous experiments on the same topic [37[ll"12|- The 
scattered intensity at a given wave vector Q is measured 
in the spin-flip (SF) and non-spin-flip (NSF) channels, 
with the neutron spin polarization H//Q. For this polar- 
ization, the full magnetic intensity always appears in the 
SF channel. Due to a large neutron depolarizion when 
entering the SC state in Bi2212 samples, the identifica- 
tion of any magnetic signal below Tc is prohibited using 
our current neutron polarization device [4] . 



Follo wing pre vious studies in bilayer Y123 [ll|-|l4| and 



Hgl201 [15| 



; prev 

, ll6| , the search for a long range magnetic or- 
der in the PG phase is performed on Bragg reflections 
(1,0,L)-(0,1,L) with integer L values. The scattered in- 
tensity in the SF channel on a Bragg reflection (Isf) is 
dominated by the leakage of the NSF intensity into the 
SF channel, whose magnitude gives the bare flipping ratio 
(FR°(T)), characterizing the neutron beam polarization 
quality and stability. On top of this signal, the intrinsic 
magnetic response (Imag)7 of much weaker intensity, can 
develop once a magnetic order settles in below a certain 
temperature. The scattered intensity in the SF channel 
then reads: 



Isf = Insf/FR°{T) + I„ 



(1) 



As previosuly discussed [1J|, the neutron intensity / 
stands in Eq. [1] for the intrinsic Bragg intensity, i.e. 
the intensity measured at the Bragg peak to which a 
background is removed. The background is typically 
measured away from the Bragg position, for instance at 
Q=(0.9,0,L) for Q=(1,0,L). 
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FIG. 2: (color online) Variation of the magnetic signal mea- 
sured above Tc at Q=(1,0,L) as a function of the wave evec- 
tor perpendicular to the Cu02 plane: Ql = for differ- 
ent cuprates: Hgl201 (UD-61) ;iH|, Y123 (UD-54) !12], Lal24 
(UD-22) 17] and Bi2212 (UD-85) (see text for scaling proce- 
dure). The upper horizontal scale indicates the full momen- 
turn Q = V(f H)2 + (fL)2. 



The averaged crystal structure of Bi2212 is usually de- 
scribed by a Bb2b space group: in our tetragonal nota- 
tions, (H,0,L) Bragg reflections are thus observable only 
for H+L even. We then first search for magnetic scatter- 
ing at Q=(1,0,L) for odd L values. Fig. [T] shows the SF 
and NSF intensities for Q=(1,0,1) where an enhancement 
of Isf occurs below ^ 230 K indicating the appearance 
of a magnetic order in our underdoped Bi2212 sample. 
The same kind of magnetic signal appears below ^230 K 
at L=l and L=3. Increasing further L to 5 and 7, the 
magnetic signal vanishes. Searches for the existence of a 
magnetic signal at even integer L values (L=0,2) or for 
non integer L have remained unsuccessful. This study 
of the magnetic scattering along (lOL) emphasizes that 
the magnetic order is likely resolution-limited, i.e. long 
range. It should be noted that Bi2212 actually exhibits 
an incommensurate crystal structure with very strong 
satellite reflections^]. That precludes the determina- 
tion of the absolute values of the magnetic moments as 
nuclear Bragg scattering intensities cannot be estimated 
with enough accuracy. 

As in monolayer Hgl201 and bilayer Y123 [111- 



141 ] , a 3D magnetic order then develops in the underdoped 
phase of bilayer Bi2212. In monolayer Lal24 [l^, this 
magnetic order remains 2D and the magnetic scattering 
spreads along the (1,0, L) rod. The magnetic intensities 
along the c axis are reported on Fig 15] for these 4 cuprates 
families. They have been all rescaled to unity for L=0 
(In Bi2212, they are actually normalized to 0.6 at L=l, 
Ql=0.20 A~^, as the (1,0,0) peak is absent). Interest- 
ingly, Iinag exhibits a similar decay with Ql for all sys- 
tems suggesting a common magnetic origin. It is worth 
to mention that Imag decays much faster than the squared 
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FIG. 3: Inverse flipping ratio, 1/FR(T), defined as Bragg 
intensity in SF channel normalized by the Bragg intensity 
in the NSF channel: Q=(1,0,1) (full squares) and Q=(2,0,0) 
(open circles) for 3 samples: a) UD-85, b) OD-87 , c) OD-70. 



magnetic form factor of Cu^"*" and ions (reported in 
Fig. [5] for a sake of comparison). 

We now focus on the hole doping dependence of the 
magnetic signal through a comparative study with 2 
other overdoped (OD) single crystals: OD-87 js^l and 
OD-70 3l|. For a comparison of the same measurements 
in samples having different masses, it is quite conve- 
nient to plot the inverse- flipping ratio 1/FR(T) where 
the magnetic signal should appear on top of the inverse 
of the bare flipping ratio, 1/FR°(T) Fig. |3] compares 
1/FR(T) in the 3 different samples at 2 different Bragg 
spots. The weak temperature dependence of 1/FR°(T) is 
determined by an extra measurement at the Bragg peak 
(2,0,0), i.e at large |Q| where the magnetic signal be- 
comes vanishingly small and can be ignored. In con- 
trast, the inverse flipping ratio at Q=(1,0,1) shows an 
enhancement upon cooling down (Fig. |3]): such an en- 
hancement highlights the lUC magnetic order, given by 

230 



mag 



[1/FR(T) - 1/FR°(T)], that develops below T 
K in sample UD-85 and below Tmag ~ 170 K in sam- 
ple OD-87. At larger hole doping, in sample OD-70, no 
magnetic signal can be detected in the normal state. It 
is worth mentionning that we observe in Bi2212 a signal 
around optimal doping that was not possible to detect 
previously in Y123 [lH or Hgl20l[l5|. The nuclear Bragg 
(1,0,1) intensity is actually weaker in Bi2212 allowing to 
observe smaller magnetic moments. 

Fig. |4] summarizes the variation of the magnetic order- 
ing temperature Tmag as a function of the hole doping 
p, given by the empirical relationship: Tc(p)/T™^'^ = 
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FIG. 4: (color online) Hole doping dependence of the super- 
conducting critical temperature Tc (Black dots) and magnetic 
ordering temperature Tmag (blue dots) corresponding to the 
large single crystals used for the polarized neutron scattering 
study. The data are compared with the temperature depen- 
dence of Tc (open squares) and the pseudo-gap temperature 
T* (red triangles) from ab-resistivity measurements on thin 
films fsl. 



1 - 82.6(p - 0.16)2, m with T™^'^ set to 91±2 K. T^ag 
decreases continuously upon increasing p and is likely to 
vanish upon approaching a critical hole level of ~0.19, the 
end point of the PG state according to electronic specific 
heat measurements [s^. Ti„ag(p) can be compared with 
T*(p) obtained from in-plane resistivity measurements 
[3^ . At this step, it is worth pointing out that in-plane 
resistivity measurements are performed on thin films and 
T™^'^ does not exceed 82 K 



34| . One should therefore 



be particularly careful when comparing large single crys- 
tals and thin films, for which crystal growth techniques 
and/or annealing conditions are different. The compari- 
son of Tniag(p) with T*(p), reported in Fig SI shows that 
they both exhibit a similar evolution as a function of hole 
doping. In particular, T*(p) does not seem to decrease 
linearly with increasing p, but rather displays a shoul- 
der close to optimal doping (p=0.16). That persistence 
well above Tc near optimal doping is typical of the phe- 
nomenology of the PG state in Bi2212. This qualitative 
comparison suggests that the observed magnetic order 
and the PG state are likely bound to each other. 

In the absence of a common theory for the pseudogap 
to guide data analysis, each experimental technique has 
developped their own definition to determine a temper- 
ature T"^. As a result, values of T* may vary from one 
technique to another. It is therefore more meaningfull 
to consider the decay rate of T* upon increasing hole 
doping. According to ARPES "sE], Knight shift NMR 
[36] and tunneling [s^j measurements, the decay rate of 
T* and the PG energy Apc is about 1.5-1.7 from the 
weakly UD regime (Tc=82± 3 K) to the weakly OD 
regime (Tc=85± 3 K). Consistently, Tmag and the or- 
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dered magnetic moment at 100 K (|M| = ^Imag) decay 
by factors ~ 1.4 and ~ 1.7, respectively (Fig. [3]). 

The observed 3D magnetic order does not appear to 
change upon structural properties of each hole-doped 
cuprate. It does not depend neither on the number of 
Cu02 planes per unit cell nor on the nature of stack- 
ing of Cu sites along the c axis. Indeed, Bi2212 has a 
body centered structure, at variance with Y123. Hgl201 
has a simple tetragonal structure, Y123 an orthorhombic 
one along the Cu-0 bonds as well as Bi2212 but along 
diagonals. In all cases, the observed order preserves the 
crystal lattice translation invariance, corresponding to an 
IUC-(antiferro)magnetic orderQ. 

The existence of a magnetic order signals shows that 
time reversal symmetry is broken. This property in the 
PG state of Bi2212 was first inferred from the observa- 
tion of dichroic effect in ARPES measurements but 
this measurement has been the subject to a long standing 
controversy (see 6.17. fsS*] and related discussions). Keep- 
ing in mind that ARPES measurements are performed 
on thin films [ij and polarized neutron measurements on 
large single crystals, the onset of time reversal break- 
ing symmetry found in both types of measurements for 
samples UD-85 is in a very good agreement (~200K for 
ARPES and ^230 K for polarized neutron diffraction). 
The doping dependence obtained in both types of exper- 
iments also match each other quite well. 

Symmetry breaking within the PG state is now cor- 
roborated on various cuprates by different techniques [Sj 
in addition to polarized neutron diffraction and the cir- 
cularly polarized ARPES measurements. In monolayer 
Pbo.55Bii.5Sri.6Lao.4Cu6+5[l|, ARPES, polar Kerr effect 
and time-resolved reflectivity measurements performed 
on the same single crystal demonstrate that a phase tran- 
sition breaking time reversal symmetry takes place at T*. 
To further characterize the actual order parameter, it is 
interesting to notice that the existence of an lUC order 
associated with the PG state in Bi2212 system is also 
supported by recent analyses 0] of scanning tunneling 
microscope (STM) images. An electronic nematic order 
has been observed as the 90° rotational symmetry is bro- 
ken in the Cu02 unit cell, yielding distinct electronic 
densities on oxygen sites along a and b crystal axis. In 
principle, this lUC order differs from the CC-loop order. 
However, the mean-field analysis of the different lUC- 
ordering possibilities in the 3-band Emery model indicate 
that the electronic nematic order and the CC-loop order 
could actually coexist [39]. 

Finally, the concomittant observations in Bi2212 of a 
dichroic effect at the antinodal wave vectors and a mag- 
netic signal on Bragg reflections (1,0, L) point towards a 
PG order parameter breaking inversion symmetry 0, [l^j 
with the result of dismissing the model of spin moments 
on oxygen atoms ll| and promoting the loop current 
CC-6'n phase 0. Recent observations of two nearly dis- 
persionless magnetic excitations in IIgl201 by polarized 



inelastic neutron scattering [40[ give actually further evi- 
dence for the presence of a such discrete order parameter 
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